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We provide a detailed description of a general procedure by which a nano/micro-mechanical resonator can 
be calibrated using its thermal motion. A brief introduction to the equations of motion for such a resonator is 
presented, followed by a detailed derivation of the corresponding power spectral density (PSD) function. The 
effective masses for a number of different resonator geometries are determined using both finite element method 
(FEM) modeling and analytical calculations. 



I. INTRODUCTION 

Nano/micro-mechanical resonators are fascinating devices, 
having microscopic dimensions and high quality factors (Qs) 
that allow them to be used in precise experiments. They are 
therefore ideal for many measurement applications, including 
biosensing [1-4], charge measurements [5], mass sensing [6- 
9], accelerometry [10], temperature sensing [11, 12], torque or 
force transduction [13-15], single spin detection [16], and ob- 
servation of qubits [17]. In addition, nano/micro-mechanical 
resonators provide an excellent platform for investigating sci- 
entific phenomena, such as vortices [18, 19], quantum zero- 
point motion [20-22], molecular structure [23], gravity [24], 
superfluidity [25], and the Casimir force [26]. 

As nano/micro-mechanical resonator systems become ca- 
pable of performing ever more precise experiments, the sen- 
sitivity of the detection methods used to observe the devices 
must increase commensurately. Current state of the art tech- 
niques used for this detection include capactive measurements 
[27], free-space optics [28-30], piezoresitivity [31], and op- 
tomechanics [ 1 3, 32-35] . These methods offer extremely sen- 
stive transduction, achieving sub-attometer displacement sen- 
sitivity [36], yoctogram mass resolution [9], and force mea- 
surements on the order of attonewtons [15]. 

In order to quantify such sensitive measurements, it is cru- 
cial that the calibration of these devices be performed ac- 
curately. As the dimensions of nano/micro-mechanical res- 
onators continue to decrease, this can prove to be a difficult 
task. However, there exists a powerful method by which 
the motion of these resonators can be calibrated via obser- 
vation of their random thermal motion. This is known as 
thermomechanical calibration and is performed by invoking 
the equipartition theorem, which establishes a relationship be- 
tween the thermal energy of a device and its mean squared 
amplitude of motion [37]. Many of the detection methods 
mentioned above are sensitive enough to measure the thermal 
motion of devices on the micro or nanoscale. In addition, it 
has been shown that this calibration scheme agrees well with 
other methods, such as physically pushing on the device [38], 
but it has the advantage of being noninvasive. 

Thermomechanical calibration has been used to measure 
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the motion of many systems, including cantilevers [29, 30, 33, 
38^3], doubly clamped beams [10, 44-46], strings [28, 32], 
torsional resonators [5, 13, 47-49], and membranes [50, 51]. 
However, to the best of our knowledge, no single standalone 
document exists that describes the method of thermomechan- 
ical calibration in its entirety. We have noticed that a num- 
ber of mutually inconsistent approaches to performing such 
a calibration can be found in the literature. This manuscript 
intends to elucidate the process of thermomechanical cali- 
bration by providing a full, self-consistent derivation of the 
procedure and ensuring that all pertinent equations are ex- 
plained and defined. In this way, a general process is de- 
veloped by which nano/micro-mechanical resonators can be 
calibrated — the only difference between varying geometries 
being reduced to a single parameter, the device's effective 
mass. In this paper, the effective mass is calculated both an- 
alytically and numerically, for a number of common device 
geometries. These calculations are then benchmarked against 
the well-known methods of Sader et al. [52] and Cleveland et 
al. [53], which have been developed to calibrate atomic force 
microscopy cantilevers. In addition, we extend the methods of 
Sader et al. and Cleveland et al. to encompass other resonator 
geometries, such as doubly clamped beams and torsional res- 
onators. 

This document — which includes both a general overview of 
the process and specific examples — provides a clear, straight- 
forward recipe for caiTying out thermomechanical calibration 
for any nano/micro mechanical resonator system, regardless of 
its complexity. 

n. MOTION OF A RESONATOR 

We begin by describing the motion of an extended resonator 
structure in the linear regime. In general, the complete, three- 
dimensional motion will be given by a displacement function 
R(x,f ) that can be broken down into an infinite number of in- 
dependent modes. To each mode we ascribe a label, n, allow- 
ing for the displacement function to be expressed as [44, 54] 

R(x,f)=l^fl„(f)r«(x), (1) 

n 

where r„(x) is the mode shape of the n* mode, and a„[t) is a 
function that describes the time dependence of the resonator's 
motion. Here we choose to normalize r„ (x) such that the max- 
imum value of |r„(x)| is unity. This choice of normalization 
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ensures that a„{t) has units of distance [as the normalization 
leaves r„(x) unitless] and is in 1-1 correspondence with the 
resonator's true physical displacement at the point of measure- 
ment. It should also be noted that by normalizing this way, the 
standard convention of orthonormality between the modes is 
lost. 

Many devices can be described by simple models, for which 
we are able to determine r„ (x) analytically. When the device's 
motion is more complicated, however, we turn to FEM mod- 
eling to determine the mode shape of the structure. In this 
paper, these two methods are compared in a number of situ- 
ations where the analytical mode shapes of the resonator can 
be determined. 

The function a„{t) is determined by mapping each of the 
resonator's independent modes to a damped harmonic oscilla- 
tor of the form [42] 

2 Fit) 
an + —a„ + (0„a„ = . (2) 

Qn n%f{.„ 

Here, o,,, Qn, and meff „ are the angular frequency, quality 
factor, and effective mass for the n* mode of the resonator 
F{t) is a time-dependent driving force. An effective spring 
constant for the rfi^ mode can also be defined and is given by 
^eff.n = 'Weff.nft',? [42]. A general expression for mgff „, along 
with its determination for a number of simple models, will be 
outlined in Sec. V. 

Often a nano/micro-mechanical resonator's motion can be 
reduced to a one-dimensional displacement function 

Zix,t) ^Y,Zn{x,t) =Y,a„{t)Un{x), (3) 

n n 

where m„(x) describes the one-dimensional mode shape for 
the n* mode of the resonator as a function of position x, and 
a„(f) is unchanged from Eq. (I). As in the three-dimensional 
case, we choose to normalize u„{x) such that the maximum 
value of |m„(x)| is unity. This equation will be used when the 
displacement of a nano/micro-mechanical resonator is simple 
enough to be modeled by a one-dimensional function. Equa- 
tions (1) and (3) fully describe the motion of these devices in 
the linear regime and will be used when explicating the ther- 
mal calibration process. 

III. POWER SPECTRAL DENSITY 

The various methods used to detect the motion of a 
nano/micro-mechanical resonator share one major similarity: 
the motion of the resonator in question is encapsulated by 
some time-varying signal, typically a voltage. Generally, this 
measured signal is transformed into the frequency domain, 
providing a spectral representation of the resonator's motion 
that is peaked (up to small corrections in \IQ) at its resonance 
frequency. By squaring this spectrum and dividing by the res- 
olution bandwidth of the instrument, the measured one-sided 
power spectral density ^-^(a)) of the signal is produced [50]. 



When thermally calibrating a mechanical resonator, we 
must be able to relate the device's time averaged motion to 
its thermal energy. To do this, the equipartition theorem is 
used to express the device's mean-square amplitude of motion 
{a^it)) in terms of its thermal noise PSD. In this section, we 
derive a theoretical expression for the PSD corresponding to a 
damped harmonic oscillator undergoing random thermal mo- 
tion. This function can be used to properly calibrate thermal 
PSD data for any nano/micro-mechanical resonator 

We introduce the autocorrelation function R^zit), which de- 
scribes how the signal a„{t') is related to itself at a later time 
t' + t and is defined as [55] 

1 /■^o 

R„(t) = lim — / dt' aJt')an(t' + t). (4) 
To^oo To Jo 

The two-sided power spectral density function P^^{co), can 
then be obtained from the Fourier transform of the autocor- 
relation function [55] 

Pzz{(0)= J dtRzz{t)e''^' . (5) 

Likewise, the autocorrelation function can be produced by 
performing the inverse Fourier transform 

Rzzit) = ^ j'^JC0Pzz{(0)e-''" . (6) 

7^2(0)) is defined such that it spans all frequencies, both pos- 
itive and negative, and will give the total power of the signal 
if integrated over this entire frequency range. For real-valued 
signals, as is the case in physical situations, Pzz{(o) is an even 
function [56]. Therefore, we introduce Szz{o)), which is de- 
fined on the frequency interval from to 0° such that it pro- 
duces the power in the signal when integrated over all positive 
frequencies. From these two definitions, we find that 5,^(0)) 
and Pzzio)) are related by 5,,(a)) = 2Pzz{o}) [56]. This rela- 
tion allows the total power in the signal to be conserved when 
each spectral density function is integrated over the frequency 
interval on which it is defined. In the rest of the derivation, 
Szzico) will be used, as it provides a theoretical spectral den- 
sity that can be fit directly to physical data. 

We define the quantity known as the mean-square ampli- 
tude of our signal as [55] 

(fl„2(f)) = ^ /\f[fl„(f)]'. (7) 

Upon inspection, we see that this quantity is equal to Eq. (4) at 
zero time shift, provided that we have sampled the signal a„ (f ) 
over a sufficiently long time window (i.e., To >> 27r/a)„). Us- 
ing Eq. (6), along with the relation between between Pzzi^'i) 
and Szz{(o), we can then express {a„^{t)) in terms of its PSD 
via [42] 

(a„2(f)) f d(0Szz{(0). (8) 

27: Jo 
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This result provides us with the relationship between the 
mean-square amplitude and PSD for a„{t). By Fourier trans- 
forming Eq. (2), we can relate S^^{co) to (to), which is the 
PSD due to the forcing function F{t) [55]: 



S,A(o)^\xico)\^SFFico). 



(9) 



The factor of proportionality comes from the frequency re- 
sponse function [10], 



frequency, it is convenient to express S^^{(o) in terms of v. 
This is done by taking co„ — 2nv„, resulting in 



kBTVn 



2nhn,if^,Qn (v2 - v2)2 + (vv„/e„) 



(15) 



Alternatively, it is sometimes useful to express S^v) in terms 
of the effective spring constant k^^ ^ [29, 30, 41, 47, 50]: 



^{an{t)} ^ 1 

^ {F{t)} meff „ (C(j2 - ©2 + i(0(On/Qn) ' 



(10) 



Here ^{a„(f)} and ^{F(f)} are the Fourier transforms of 
a„(f) and F{t). For thermomechanical calibration, F{t) is 
generated by random thermal noise, which is broadband in 
the region of interest, producing a Sff{co) that will be flat 
with respect to frequency. We can therefore replace Sff{(o) 
by a constant thermal force noise 5*^. This factor can be de- 
termined by combining Eqs. (8)-(10), which results in 



cth 



dco 



2?r<ff,„ Jo (0)2 - (o^y + {coco„/Q„ 



(11) 



Computing this integral by contour integration, gives 
(fl„2(f)) = e„/4a)>eff,„2 [54]. 

Invoking the equipartion theorem, {a^{t)) can be related to 
the thermal energy of the resonator. The equipartition theorem 
states that every degree of freedom that contributes a quadratic 
term to the total energy of a system has an average energy of 
ksT /2 [37]. For the time averaged potential energy of the one- 
dimensional harmonic oscillator considered here, this theorem 
is mathematically represented as [39] 



(12) 



It can be seen by investigation of Eq. (1) that the magnitude of 
a„(f) is arbitrarily determined by our choice of normalization 
for r„(x). Therefore, our effective mass must be chosen very 
carefully to match our definition of a„(f) in order to preserve 
the validity of the equipartition theorem. This subtlety will be 
discussed in detail in Sec. V. 

Combining the result from Eq. (11) with the equipartition 
theorem gives [42] 



ctn 

^>ff — 



Qn 



(13) 



allowing us to express the oscillator's PSD as [10, 42, 57] 



4kBT(X)„ 



meff.nQn 



(14) 



However, since the frequency we read off our instrument is 
given in terms of a natural frequency, v, and not an angular 



2kBTvl 



^keff^„Q„ 



(V2 



(16) 



Equations (14)-(16) provide the power spectrum for the ther- 
mal noise of our resonator However, we expect a number 
of other sources of noise to be present in our detection sys- 
tem, such as detector dark current and optical shot noise [45]. 
These signals will combine with the thermal noise of the res- 
onator to produce the power spectrum read out for the un- 
driven device. If we assume these other sources of noise to be 
white, we can express our total PSD as [50] 



Svv{v) = S'^y + aS,.{v) 



(17) 



where Syy is a constant offset to account for the white noise 
of the detection apparatus and o: is a factor with units V^/m^ 
to allow for conversion of our theoretical displacement PSD 
(in units of m2/Hz) to the measured voltage PSD (in units of 
V2/Hz). In other words, a tells us how efficiently the detec- 
tion system converts the motion of the resonator into the volt- 
age signal that is measured and is closely tied to the sensitivity 
of the detection system. 

The experimentally measured PSD is given as 5^"^ (v) = 
y (v)^/Av [28] where V{v) is the voltage spectrum of the res- 
onator's motion obtained with a measurement bandwidth Av. 
The experimental PSD can be fit with Eq. (17) near a thermal 
resonance using Syy, a, V„, and Qn as fit parameters. Once a 
has been determined, it can be used to convert the PSD data 
S'^y from units of V2/Hz to m2/Hz, the result of which we 
will call ^x?. The resonator's displacement spectrum is then 
determined from \/^^, given in units of m/ VHz [29]. From 
this data, we can determine the maximum displacement of the 
resonator, as well as the displacement sensitivity of the detec- 
tion method, the latter of which is given by ^yS^y/a. The 
displacement sensitivity is often used as a figure of merit for 
a detection system [29, 32, 33], as it is a measure of the min- 
imun detectable signal for that setup. 

It should be noted that a is highly dependent on a num- 
ber of factors, for instance photodiode gain or probe position, 
and therefore an a value extracted from a specific set of data 
should be associated only with that data set itself and not the 
detection system. This being said, it is possible to conceive of 
an experiment where all parameters are carefully controlled 
such that the a value for repeated data sets would be simi- 
lar and could be combined to determine an average a value. 
Similarly, a number of data sets originating from a carefully 
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controlled system could be combined into an averaged sig- 
nal and fit with Eq. (17) to provide a statistically significant 
value of a. The theoretical PSD functions obtained in this 
section are a general result for a damped harmonic oscillator 
excited by random thermal fluctuations and can be used for 
any nano/micro-mechanical resonator. The only difference in 
Eq. (15) [Eq. (16)] that exists between varying resonators is 
the effective mass meff„ (spring constant kgff „) associated with 
the resonator's geometry. Therefore, once we have the PSD 
function given by Eqs. (15) and (16), we need only focus on 
calculation of the device's effective mass. For this reason, the 
calculation of the effective mass for a number of resonator ge- 
ometries will be the focus of the remaining portion of the doc- 
ument. For simplicity, we deal only with the quantity Weff,,,, 
as ^eff.n is directly related to /«eff,„ through the mode's experi- 
mentally determined resonance frequency. 



where L is the length of the beam and A„ = L (^pAco^/El) ^''^ 
is a dimensionless parameter This differential equation has 
the general solution [58] 



Unix) =:A„sin ( +B„cos f ^-"^ 



C„ sinh ( 



D„ cosh ( 



(21) 



Enforcing the beam's boundary conditions allows us to deter- 
mine the values of X„, as well as A„, B„, C„, and D„ (up to 
a normalization constant), thus solving for the mode shape of 
each of the beam's n modes. It is then possible to calculate the 
frequency corresponding to the mode of the beam using 



IV. ANALYTICAL THEORIES OF SIMPLE RESONATORS 

In general, the effective mass associated with the specific 
mode of a resonator will depend on its mode shape. For this 
reason, it is useful to investigate analytical expressions for the 
mode shapes and resonance frequencies of some simple ge- 
ometries that are commonly used for nano/micro-mechanical 
resonators. 



A. Analytical theory of beams 

A large number of nano/micro-mechanical resonators pos- 
sess beam-like geometries and can be effectively modeled as 
macroscopic beam-like structures. The general equation of 
motion for a long thin beam of uniform cross section that is 
oscillating in one dimension is [58] 



EI-. 



pA 



--F{x), 



(18) 



where E, I, p, and A are the Young's modulus, moment of 
inertia, density, and cross-sectional area of the beam. F{x) is 
an arbitrary position-dependent loading function. We know 
that z{x,t) will be in the form of Eq. (3), so we can solve 
Eq. ( 1 8) for each of the resonator's n modes separately. The 
oscillatory motion of the beam is encapsulated by fl„(f) and, 
for high Q resonators, is given approximately as 



2n InL'^ V pA ' 



(22) 



A method of determining a cantilever's effective mass by 
applying a load to it and observing its deflection has been de- 
veloped by Sader et al. [52]. It is therefore useful to determine 
the deformed mode shape w{x) when a beam is subjected to a 
static loading function F{x), which can be found using [58] 



(fiu , , 

EI^_=-M[xl 



(23) 



where M{x) is the bending moment of the structure due to 
the specified loading F{x). This can be solved analytically 
for a number of different loadings, some of which are very 
important in defining effective masses and will be discussed 
below. 

Solving the above equations for the specific cases of a 
singly clamped beam (cantilever) and a doubly clamped beam 
will be very useful in providing an analytical value for the 
effective mass in these two simple cases. In the subsequent 
sections, we will treat each of these cases separately, deter- 
mining both dynamic and static displacements. These will be 
used later to determine the effective mass of each device. 



1. Cantilevers 

The singly clamped beam, or cantilever, is subject to the 
boundary conditions [58] 



a„{t) = c„cos(a)„f) +ii„ sin(a)„f), (19) 

where c„ and d„ are mode dependent constants. Inputting 
z„(x,f) into Eq. (18) with this choice of fl„(f), we obtain for 
the free motion {F{x) = 0] of a beam, a differential equation 
for the mode shape of the mode 



These boundary conditions tell us that the cantilever is fixed 
and flat at x = 0, while free at jc = L, as can be seen in Fig. 1 . 
Inputting these conditions into Eq. (21), we find that X„ obeys 
the relation [59] 

cos A„ cosh A,„ + 1=0, (25) 
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Position Along Device (x) 

FIG. 1: (a)-(f) Mode shapes for the first six modes of a cantilever. 
The solid (black) and dashed (green) lines respresent the positive 
and negative maximum displacement (i.e., 180° out of phase in time 
with respect to each other). Each mode is normalized such that its 
maximum displacement is equal to one. 



the solutions of which are summarized in Table 1. The mode 
shape of a cantilever is then determined to be [54] 



k = 



3EI 

Wmax 



(29) 



which will be used later to calculate the effective mass of a 
simple cantilever using the method of Sader et al. 

2. Doubly clamped beams 

The boundary conditions for a doubly clamped beam are 
[58] 



M„(0) = u„{L) 



du„ 
dx 



(0) 



du„ 
dx 



(L)=0. 



(30) 



Physically, this means that the beam is both flat and fixed at 
each of its endpoints, as can be seen in Fig. 2. Inputting these 
conditions into Eq. (21), we obtain the following condition on 
Xn for doubly clamped beams [59] 



cos A„ cosh A„ — 1=0. 



(31) 



The solutions to this equation are summarized in Table I. The 
mode shape for a beam can also be solved and is found to have 
the form [59] 



Un{x) = 



Cn 



cosh 



— cos 



A;^J 



sinh 



XnX 

1^ 



— Sin 



(26) 



where K„ = 2(sinA„coshA„ — cos Ai sinh A„), C„ = sinhA„ + 
sin A,,, and Sn — coshA„ +cosA„. Notice that our choice of 
normalization constant ensures that u„{L) — 1. This satisfies 
the mode shape normalization condition found in Sec. II, as 
the maximum amplitude of the cantilever is always at its free 
end. The first six normalized mode shapes of a cantilever can 
be seen in Fig. 1 . 

It is also useful to determine the deformed mode shape w{x) 
of a cantilever when a static force Fq is applied perpendicular 
to its free end [F{x) = Fo5{x — L)], so that its effective mass 
can be approximated using the method of Sader et al. By cal- 
culating the bending moment M{x) for such a loading force, 
we can solve for w(x) using Eq. (23), resulting in [59] 



w{x) 



Fox' 
2EI 



(27) 



The maximum displacement of the cantilever due to this load- 
ing will clearly occur atx = L, so that 



w{L) 



3EI 



(28) 



This allows us to define a corresponding spring constant 



Unix) 



K„ 



cosh 



X„x 

1^ 



-s„ 



sinh 



XnX 



— Sin 



X/^x 



, (32) 



where K„ is a normalization constant chosen such that the 
maximum value of |m„(x)| is 1, C„ = sinhA,, — sinA,,, and 
S„ ~ coshA,, — cosA„. It is interesting to note the similari- 
ties between Eqs. (31) and (32) for a doubly clamped beam 
and the corresponding Eqs. (25) and (26) for a cantilever 

Unfortunately, the position at which |M„(.5t)| is maximized 
is not obvious in the case of the doubly clamped beam, aside 
from the fundamental mode where the maximum is located 
at X = L/2. For this reason, the determination of K„ must be 
carried out numerically. The first six normalized mode shapes 
for a simple doubly clamped beam can be seen in Fig. 2. 

As with the cantilever, we are also interested in the static 
mode shapes of doubly clamped beams due to specific load- 
ings F{x). This allows us to extend the method of Sader et 
al. to approximate the effective mass of this structure. For a 
point loading force Fq applied perpendicularly to the center of 
the beam [F{x) = Fq8{x — L/2)], Eq. (23) can be solved to 
determine the deflected mode shape w{x) as [59] 



w{x) 



Fqx / L X 
4EI V4 ^ 3 



(33) 



The maximum displacement occurs atx = L/2 and is found to 
be 
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Position Along Device (x) 



FIG. 2: (a)-(f) Mode shapes for the first six modes of a doubly 
clamped beam. The solid (black) and dashed (green) lines respre- 
sent the positive and negative maximum displacement. Each mode is 
normalized such that its maximum displacement is equal to one. 




L L 

Position Along Device {x) 



FIG. 3: (a)-(f) Mode shapes for the first six modes of a string. The 
solid (black) and dashed (green) lines respresent the positive and neg- 
ative maximum displacement. Each mode is normalized such that its 
maximum displacement is equal to one. 



Mode A„ 

Number (n) Cantilever Doubly Clamped Beam 

1 TSTSl 4.7300 

2 4.6941 7.8532 

3 7.8548 10.9955 
n>3 {n- l)7i: + n/2 nn + n/l 



TABLE I: Values of A„ for cantilevers and doubly clamped beams. 
The last row provides an approximation for A„ for n> 3, which can 
be obtained by the fact that for large n, coshA„ is so large that we 
require cosA„ to be very close to zero to satisfy Eqs. (25) and (31). 



modify Eq. (18) through the existence of a stress dependent 
term giving [58] 

d'^Z d^z d^z 

where a is the axial stress in the beam with the axial tension 
in the beam due to this stress given by = aA. In the limit of 
very large tension in the beam, we can neglect the first term in 
Eq. (36) and by setting F{x) = 0, we get the equation for the 
free motion of a string 



Wmax = w{L/2) 



FqL' 
192EI' 



(34) 



The corresponding spring constant k due to this force is given 
by 



Fo 192EI 



L3 



(35) 



0. 



(37) 



This equation is solved by inputting Eq. (3) for z{x,t) and 
solving for the mode shape of each of the n independent 
modes, which are described by 



0. 



(38) 



Fixing the string at its endpoints m„(0) = u„{L) ~ [28], we 
get solutions [6] 



B. Analytical theory of strings 

Another device of great interest is the nanostring [6, 11, 28, 
60]. The advantages of such structures are that they possess 
very high quality factors and can be modeled as macroscopic 
strings, which possess a very simple mode shape [28]. The 
physical realization of a nanostring is to apply tension to a 
doubly clamped beam, which can be performed by fabricating 
the beam out of materials with a high intrinsic stress, such 
as Si3N4 [60]. Applying an axial stress in this manner will 



u„{x) 



(39) 



/n7tx\ 

where |m„(j[:)| conveniently already has a maximum value of 
unity due to the nature of sin(x), satisfying our normalization 
condition. The mode shapes for the first six modes of a string 
are shown in Fig. 3. The resonance frequencies of a simple 
string can also be determined and are given as [28] 



(40) 



C. Analytical theory of membranes 

A membrane is one of the few two-dimensional 
nano/micro-mechanical systems that can be solved analyti- 
cally. The physical realization of a membrane can be thought 
of as the two-dimensional analog of the nanostring men- 
tioned above, in which a very thin plate is fabricated from 
a high stress material, such as Si3N4 [51]. Nano/micro- 
mechanical membranes are used in membrane-in-the-middle 
optomechanics experiments [51, 61, 62] and have also been 
used as electrometers with charge sensitivities below the sin- 
gle electron level [50]. In this section, we introduce the mode 
shapes and resonance frequencies of rectangular and circular 
membranes, which will allow us to compute their effective 
masses analytically in the following section. 
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1. Rectangular membranes 

The equation of motion for a rectangular membrane can be 
obtained by generalizing the equation of motion for a string, 
given by Eq. (37), to two dimensions, resulting in [63] 

To solve this equation, we break z{x,y,t) into an infinite num- 
ber of m and n modes in the x and y dimensions, respectively, 
giving a solution of the form 

z{x,y,t)=Y^Y^a„,„{t)^l,nn{x,y), (42) 

m n 

where (-^jy) = u,„{x)vn{y) is a two-dimensional separa- 
ble mode shape function that is normalized such that its max- 
imum value is unity. Upon inspection it can be seen that 
this solution is the two-dimensional version of Eq. (1), with 
r„(x) = '^mn{x,y)- Similar to the one-dimensional analysis, 
cimn{t) describes the oscillatory motion of the membrane and 
can be expressed as 



FIG. 4: (a)-(h) Mode shapes for the first eight modes of a rectangular 
membrane. Each mode is labelled with its mode index (m,n). (i) 
Diagram indicating the coordinates and lengths of each side of the 
rectangle and Ly. 

This results in the solutions u,„{x) — sm{7Zmx / L^) and v„{y) — 
sm{7Zny / Ly) , which produce the following expression for the 
mode shape of the membrane [63]: 



Vmn (x, y ) = sin 1^ -j— Y~ j ■ ^^^^ 

Conveniently, the solution in terms of sin functions gives 
Wmn{x,y) the normalization condition that its maximum value 
be unity. The first eight mode shapes for a rectangular mem- 
brane are shown in Fig. 4(a)-(h). Using this result, we can find 
the resonance frequency of each of the (m,n) modes [63]: 




amr,{t) = c,„„cos(a),„„f) +£/„,„ sin(a),„„f) 



(43) 



2. Circular membranes 



Inputting this function for a„m{t) into Eq. (41), we obtain an 
equation of motion for the mode shape of the membrane 



~dx^ 



' dy^ 



0. 



(44) 



We can solve this equation using the boundary conditions for a 
rectangular membrane of length along the jc-axis and length 
Ly along the y-axis [see Fig. 4(i)] fixed at all of its edges [64], 
which are given by 



Wmn{0,y) = Wmn{x,Q) = Wmn{Lx,y) = Wmn{x,Ly) = 0. (45) 



Circular membranes, used in a number of experiments [65- 
67], are another interesting nano/micro-mechanical resonator 
that have an analytically solvable mode shape. 

We begin this analysis by first noting that Eq. (4 1 ) can be 
generalized to a coordinate independent form by introducing 
the two-dimensional Laplacian operator, V^, which gives [63] 



(48) 



We note that for Cartesian coordinates, = /dx^ + 
d^/dy^. The equation of motion in this form can now be 
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solved using any system of coordinates we prefer, provided 
we apply an appropriate definition for V^. 

When working with circular membranes, it is convenient to 
introduce cylindrical coordinates, as it mirrors the symmetries 
of the system. Transforming Eq. (48) to cylindrical coordi- 
nates using the appropriate Laplacian operator gives [63] 



ld_ 
s ds 



dz 



1 dh 



-0, 



(49) 



where we have introduced the radial coordinate s and the an- 
gular coordinate as defined in Fig. 5(i). We solve this equa- 
tion of motion by inputting a solution of the form 



(50) 



Here \i/,„„{s,(j)) — S{s)<P{(j)) is given as the product of the func- 
tions S{s) and 't>(0), describing the membranes motion as a 
function of s and (j). We can then input Eq. (50) into Eq. (49), 
producing the equation of motion for the mode shape of the 
membrane 



1 d / dS 



1 d^<P 



50 = 0. 



(51) 



This equation can be solved for a fixed circular membrane of 
radius a by ensuring that the boundary condition \i/,n„{a, 0) — 
is met [64]. As well, we require the periodic condition that 
Wmn{s,(j>) = y^mn{s, (j> +2n7z), which ensures that y/mni-^, 0) is 
single-valued. The solutions resulting from these conditions 
are 't>(0) = cos(m0) and S{s) = Jm{a,„„s/a), where is the 
Bessel function of the first kind and a,„„ is determined from 
the solutions of /,„ («,„„) = 0. Table 11 contains the first few 
values of «,„„. Note that the general solution of S{s) would 
also contain a Bessel function of the second kind, however, 
this function diverges as 5 ^ 0, so it is discarded for physical 
reasons. As well, we note that 4>(0) = sin(m0) is also a pos- 
sible solution to Eq. (51), corresponding to a n:/2 phase shift 
in 't>(0) meaning each mode with m > is two-fold degener- 
ate. However, both of these degenerate modes have identical 
effective masses, so we concern ourselves only with the so- 
lution 't>(0) = cos(m0) in this document. The mode shape 
function for a circular membrane is then given by [63] 

Wmn {s,<P)= K„, COS (ot0 )7,„ ^ , (52) 

where K,n is a normalization constant chosen such that the 
maximum value of \j/,n„{s,(l)) is unity. The first eight modes of 
a circular membrane are shown in Fig. 5(a)-(h). We note that 
due to the properties of the Bessel functions, radially sym- 
metric modes (m = 0) have their maximum value at s — 0, 
while all other modes (m > 0) have a node at 5 = [see 
Fig. 5(a)-(h)]. Therefore, we have that Kq = \, which is a 
consequence of the fact that yo(0) = 1. For m > 0, K,n is given 
by [l/Zmlimax)] ■, where imax is the radial position at which 







n = 1 


n = 2 


n = 3 


;i = 4 


n = 5 




— 


9 4.04.0 


5.5201 




1 i . / 7 1 J 


14 Q^flQ 


m 


= 1 


3.8317 


7.0156 


10.1735 


13.3237 


16.4706 


m 


= 2 


5.1356 


8.4172 


11.6198 


14.7960 


17.9598 


m 


= 3 


6.3802 


9.7610 


13.0152 


16.2235 


19.4094 


m 


= 4 


7.5883 


11.0647 


14.3725 


17.6160 


20.8269 


m 


= 5 


8.7715 


12.3386 


15.7002 


18.9801 


22.2178 



TABLE II: Values of the n''' zeroes for Bessel functions of the ra''' 
kind, OCmn, for m = — 5 and n = 1 — 5. These values are taken from 
[68]. 
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FIG. 5: (a)-(h) Mode shapes for the first eight modes of a circular 
membrane. Each mode is labelled with its mode index (m,n). (i) 
Diagram indicating the coordinates and radius a of the circle. 



Jm(s) has its maximum value. To determine this value, we use 
the relation [63] 



dJm[s) 

ds 



[Jm-\{s)-Jm+\[s)\^Q, m>0, 



(53) 



which allows us to determine the extrema of J,„ (s) via the con- 
dition J,n-i{s) = Jm+i (■s)- It turns out that the first solution of 
this equation (i.e., the one closest to zero), produces imax- Us- 
ing this relation, the first few values of K,n are determined and 
presented in Table 111. 

Finally, we find that the frequency for each of the modes 



K4 
K5 



1.0000 
2.9536 
4.2251 
5.2995 
6.2609 
7.1456 



TABLE III: Values of the normalization constant K,„ 
from Eq. (5), for m = — 5. 



[1/ Jmi^max)] 
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will be given by 

v,™ = — = — J . (54) 

In 271 \ p a 

We are now equipped with the necessary solutions required 
to analytically determine the effective masses of these simple 
structures. Such calculations will be investigated in the next 
section. 



V. EFFECTIVE MASS 



determined. It should be mentioned that since the PSD cal- 
culated in Sec. Ill was for the time-dependent displacement 
function a„{t), which corresponds to the motion of the device 
at its point of maximum amplitude by virtue of the normal- 
ization condition for r„(x), Eq. (58) corresponds to the effec- 
tive mass for a point-like measurement taken at this position 
of maximum displacement. For a point-like measurement at 
some other position on the resonator, xq, the effective mass 
defined in Eq. (58) must be modified by dividing by |r„(xo)p. 
Therefore, the effective mass is in fact a function of the posi- 
tion at which we measure the device, leading to a more general 
definition 



In order to precisely calibrate a nano/micro-mechanical res- 
onator thermomechanically, one must be careful when deter- 
mining its effective mass, as numerous different methods and 
definitions used for calculating effective mass appear in the 
Uterature [9, 32, 34, 42, 44, 52-54, 57, 69-77]. What must be 
remembered is that we are free to use any definition of effec- 
tive mass to calibrate our devices, provided that an appropri- 
ate modification is made to the measured amplitude a„ (f ) such 
that the equipartition theorem holds. In other words, in order 
to describe our resonators as damped harmonic oscillators, we 
must include a correction term in either the effective mass or 
the measured amplitude of the device (or both) that accounts 
for the extended nature of the device. However, it is impera- 
tive that one applies this correction factor self-consistently to 
ensure that the equipartition theorem holds. In this paper, we 
have chosen to define our effective mass in such a way as to 
leave the measured amplitude of motion uncorrected, that is 
to say a„(f) tracks the device's true physical motion at the po- 
sition of maximum displacement. When mentioning effective 
mass in this paper we are referring to this definition. 

The effective mass of the n* mode of a resonator can be de- 
termined by examining its potential energy. For the mode, 
each small element dV of the resonator will have the potential 
energy dU corresponding to that of a harmonic oscillator with 
a mass element dm — p {r)dV given by 



dU=^p{r)(0^\a„{t)r„{x)\^dV. 



(55) 



The total potential energy of the mode is then given by inter- 
grating over the entire volume of the device, V, to get 



U - -C0f,\an{t)\' JdVpir)\r„{x) 

1 2 2 



where we have defined [34, 69-71, 74] 



(56) 
(57) 

(58) 



With this definition, m^ff ,, can be calculated once the normal- 
ized mode shape r„(x) for the mode in question has been 



mes,,, = dVp{r) |r„(x) 



'«eff,«(xo) = 



/^Vp(r)|r„(x)p 
|r„(xo)|2 



(59) 



Due to the fact that our normalization condition requires 
|rn(xo)| < 1, we realize that measuring the device's motion 
at some point xq other than the point of maximum displace- 
ment results in an increase in effective mass. By inspection 
of Eqs. (14)-(16), it can be seen that a larger effective mass 
results in a smaller PSD signal. It is therefore advantageous 
to observe a device at its position of maximum displacement, 
in order to increase the signal to noise ratio of the measure- 
ment. For this reason, in the remainder of the document it will 
be assumed that the measurement is being performed at the 
position of maximum displacement. 

The definition given by Eq. (58) includes devices that have 
a position-dependent density, such as multi-layer structures 
[78] and photonic crystal devices [10, 22]. However, most 
resonators have a uniform density, so we can take p (r) = p in 
Eq. (58), resulting in an effective mass given by 



?Meff,« =p dV 



(60) 



This result can be further simplified in the case of one- 
dimensional motion of a device with uniform cross-sectional 
area A, in which case the effective mass will be defined as [44] 



: pA / dx |m„(ji:)| 
Jo 



(61) 



In this document, we collectively call Eqs. (58)-(61) the ef- 
fective mass integral (EMI). The EMI provides a general pro- 
cedure for determining the effective mass of any nano/micro- 
mechanical resonator As such, it is versatile and can be ap- 
plied in a number of complex situations. For instance, there 
often exist other entities, such as supports or under-etched re- 
gions, which result from fabrication processes and move with 
the device itself. Such motion contributes to the energy of the 
mode and therefore it must be included in the effective mass 
calculation. A simple model describing only the motion of 
the device itself would therefore be an incorrect assessment 
of the effective mass. However, by using FEM simulation, it 
is possible to perform the EMI which includes the motion of 
all involved entities. 
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The EMI is also useful for structures that oscillate in more 
than one dimension. For most systems that can be described 
by simple analytical models, the motion is constrained to one 
dimension. However, many systems, such as membranes [50, 
51] and modified atomic force cantilevers [79], oscillate in 
more than one dimension. Using FEM simulation, the EMI 
provides a method to determine the effective mass attributed 
to the complex three-dimensional motion of such devices. 

While the EMI is exact in the context of this paper, other 
methods for determining the effective mass for the fundamen- 
tal mode of a cantilever have also been developed [52, 53]. 
One popular method introduced by Sader et al. is used to de- 
termine the effective mass of the first mode of a cantilever in- 
directly by calculating its spring constant k [52]. This can be 
done by applying a point force to the free end of a cantilever, 
measuring its displacement and determining k via Hooke's 
Law, F — —kx. Note that this spring constant is identical to 
the one defined in Eq. (29). This analysis can be performed 
in many ways, either analytically [58], by simulation [52], or 
experimentally [38]. Once k has been determined, Weff j can 
be calculated by the relationship k — WeffjOj^. 

Another method of calculating the effective mass for the 
fundamental mode of a cantilever was developed by Cleveland 
et al. [53]. In this method, the effective mass is determined by 
adding a point mass to the free end of the cantilever, resulting 
in a reduced resonance frequency 



Vi 



(0]_ 

2k 



keffd 



(62) 



where M is the added mass. 

The two alternative methods to determine effective mass 
described above are developed for cantilevers. However, here 
we adapt these methods to use them to determine the approx- 
imate effective mass for the fundamental modes of other sim- 
ple nano/micro-mechanical resonators. 

It should be noted that these methods quickly break down 
as we move beyond the first mode or to more complex de- 
vices. It is therefore advantageous to use the EMI whenever 
possible when determining the effective masses of devices, as 
it is completely general. It should also be noted that for the 
definition of effective mass used in this paper, the EMI is ex- 
act and the methods of Sader et al. and Cleveland et al. are 
elegant approximations to this definition of effective mass. 

In the subsequent sections, we will analyze the effective 
mass ratio, defined as the device's effective mass Weff,,, di- 
vided by its geometric mass m, for the structures discussed in 
Sec. IV. The effective mass ratio is a useful, intrinsic quantity 
of a system as it is independent of size or material. 

In the following calculations that involve beam-like struc- 
tures, we have chosen to analyze geometries with a uniform 
rectangular cross section. It is important, however, to note that 
the results are identical for a beam with a thin, uniform cross 
section of any shape, provided that the structure in question 
can be appropriately modeled as beam-like. Justification for 
this claim is given by the fact that Eq. (21) was obtained for 
a beam-like structure without making any assumptions other 
than the fact that the cross section was uniform and thin. 



Mode 




Label (n) 


EMI - FEM 


1 


0.2498 


2 


0.2498 


3 


0.2499 


4 


0.2498 


5 


0.2498 


6 


0.2498 


7 


0.2497 


8 


0.2496 


9 


0.2496 


10 


0.2495 



TABLE IV: Values of meff,«/m calculated using FEM to solve the 
EMI for the first ten modes, n, of a cantilever with a uniform rectan- 
gular cross section. 



Using the methods described in this section, we calculate 
the effective mass ratio for each these structures, both analyt- 
ically and through simulation. We are then able to compare 
the values obtained using each of the methods, allowing us to 
evaluate their validity. 



A. Effective mass of a cantilever 

We begin our calculations by computing the effective mass 
ratio of a cantilever analytically by inputting Eq. (26) into the 
EMI, giving 



^jfdx \u„{x)\^ = J. 
m L Jo 4 



(63) 



This result is deceptively simple and tells us that the ratio 
m^ff ji/m is mode independent; that is, the effective mass of 
a cantilever is constant and always equal to 1 /4 of the geo- 
metric mass, regardless of mode. It should be noted that this 
result agrees with [54], provided that the mode shape func- 
tions are normalized according to the convention presented in 
this paper. 

We can also use the EMI to calculate a cantilever's effective 
mass ratio using FEM simulation. The results for the first ten 
modes of a cantilever with rectangular cross section are given 
in Table IV. It can be seen that these values agree well with 
the analytical value of 1 /4 to at least 3 decimal places. 

We continue calculating the effective mass ratio of the fun- 
damental mode of the cantilever using the methods of Sader et 
al. and Cleveland et al. We perform an analytical calculation 
using the method of Sader et al. by combining Eqs. (22) and 
(29), which can be rearranged to the form Oi = y/k/m^ff i, 
thus identifying an effective mass. The result is 



0)1 



pAL* 



IkX^ 
3 m 



meff,l 



(64) 



where we have taken ;«eff,i — 3m/Xi as our effective mass. 
Using the value of Ai found in Table 1 for a cantilever, we find 
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Method 


meff,i/m 


EMI - Analytical 


0.2500 


EMI - FEM 


0.2498 


Sader et al. - Analytical 


0.2427 


Sader et al. - FEM 


0.2425 


Cleveland et al. - FEM 


0.2438 



TABLE V: Values of mf.f{ „/m for the fundamental mode of a simple 
cantilever calculated using the different methods detailed in Sec. V. 



'Weff.i/'w = 0.2427 for this method. It should be reiterated that 
this method of determining the effective mass of a cantilever 
is limited to the fundamental mode. This can easily be seen by 
the fact that if we tried to extend the definition to the mode, 
the result would be niefi.n/ni — which incorrectly pro- 

duces a rapidly decreasing effective mass ratio as n increases. 
This is in direct contradiction to Eq. (63), which indicates that 
the effective mass ratio for a cantilever is constant. The reason 
why this method cannot be extended beyond the fundamental 
mode is that the spring constant k we calculate here is only 
meaningful for the first mode, i.e., the spring constant is mode 
dependent with our definition of effective mass. 

The effective mass ratio can also be calculated using FEM 
simulation for this method. A force is applied to the free end 
of a cantilever and, using stationary analysis, k is calculated. 
An eigenfrequency solver is then used to obtain the frequency 
of the fundamental mode of the same cantilever The effective 
mass is then determined from mgff.i — kj {7.%V\f- . Using this 
method we find that njeffj jm = 0.2425. 

Finally, we can calculate the effective mass ratio for a can- 
tilever using the method of Cleveland et al. This is performed 
using FEM simulation, in which a point mass is added to the 
end of the cantilever and the coiTesponding change in reso- 
nance frequency is observed. This analysis is performed with 
additional masses of varying magnitudes, producing the data 
points that are plotted in Fig. 6(a). This data is then fit with 
Eq. (62), from which we can extract the effective mass ratio. 
For the cantilever data, we find that mgff i/m = 0.2438 using 
this method. 

Again, this method cannot be extended past the fundamen- 
tal mode. The reason in this case is that the added mass at 
the free end, while not adversely affecting the first mode, sig- 
nificantly modifies the mode shape of higher order modes, 
preventing an accurate determination of the effective mass in 
these cases. 

It can be seen that the methods of Sader et al. and Cleve- 
land et al. provide a very good approximation to the effective 
mass for the fundamental mode of a simple cantilever. This is 
made evident by comparing mgff, i / m for the varying methods 
presented here, which are summarized in Table V. 

It should be mentioned that the effective mass calculations 
performed in this section were for a simple cantilever with 
uniform rectangular cross section, allowing us to solve for its 
effective mass analytically. In general, however, the geome- 
try of a cantilever will be more complex, requiring the use of 
FEM simulation to calculate its effective mass. 
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FIG. 6: (a) Plot of normalized frequency versus normalized added 
mass for the fundamental modes of a cantilevers (red circles), dou- 
bly clamped beams (blue squares) and strings (green triangles). The 
frequency values are normalized by dividing by the fundamental fre- 
quency of each device with M = 0. The added mass is normalized by 
dividing by the total geometric mass m of the device. Solid lines 
are fits to Eq. (62). Effective mass ratios extracted from the fits 
are meff,i/m = 0.2438, 0.3864, 0.4265 for the cantilever, doubly 
clamped beam and string, which can be found in Tables V, VII, and 
IX, respectively, (b) FEM mode shape for the fundamental mode of 
a cantilever for both the unloaded (M = 0: solid blue) and maximally 
loaded (M = m: dashed orange) cases. Similar plots for a (c) beam 
and (d) string, are also presented. 



B. Effective mass of a doubly clamped beam 

Just as with the cantilever, we will begin the effective mass 
ratio calculations for a doubly clamped beam by computing 
the analytical result obtained using Eqs. (32) and (61): 



m 



I [''dx\u„{x)\\ (65) 
L Jo 



Unfortunately, the ratio meff„/m is not mode independent for 
a doubly clamped beam as it is for a cantilever However, it ap- 
pears that for n > 4 it saturates to a value of Weff.n/'w ~ 0.4372. 
We can also compute the effective mass ratio by solving the 
EMI using FEM simulation as we did with the cantilever Both 
the analytical and FEM results are given for the first ten modes 
of a simple doubly clamped beam in Table VI. 

The method of Sader et al. was originally developed for a 
cantilever However, we can imagine an analogous method 
for doubly clamped beams in which we apply a point force in 
the middle of the beam instead of at the end, as that is where 
we expect to have the maximum displacement for the funda- 
mental mode. We can find an analytical value for mgff i/m by 
combining Eqs. (22) and (35), which results in 



\92m 



meffj 



(66) 
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Mode 
Label (n) 






EMI - Analytical 


EMI - FEM 


1 


0.3965 


0.3959 


2 


0.4390 


0.4381 


3 


0.4371 


0.4358 


4 


0.4372 


0.4353 


5 


0.4372 


0.4347 


6 


0.4372 


0.4339 


7 


0.4372 


0.4330 


8 


0.4372 


0.4320 


9 


0.4372 


0.4309 


10 


0.4372 


0.4296 



TABLE VI: Values of nieif.n/'" calculated by solving the EMI both 
analytically and using FEM for the first ten modes, n, of a doubly 
clamped beam with uniform rectangular cross section. 



Method 


raeff,i/m 


EMI - Analytical 


0.3965 


EMI - FEM 


0.3959 


Sader et al. - Analytical 


0.3836 


Sader et al. - FEM 


0.3830 


Cleveland et al. - FEM 


0.3864 



TABLE VII: Values of m^fi „/m for the fundamental mode of a sim- 
ple doubly clamped beam calculated using the different methods de- 
tailed in Sec. V. 



Here we identify ;«eff,i = I92m/X\^ as the effective mass of 
the doubly clamped beam in this equation. Using the value 
of Al for a doubly clamped beam from Table 1, we find that 
mgff i/m = 0.3836. We can also perform a similar analysis 
using FEM, which produces m^ffi/m = 0.3830. Both of these 
results can be found in Table Vll. 

We can also extend the method of Cleveland et al. for de- 
termining the effective mass of a cantilever to an analogous 
process for a doubly clamped beam. Again, instead of apply- 
ing the mass loading to the free end as with the cantilever, the 
mass is added to the center of the beam. Such a loading per- 
turbs the fundamental mode of the beam by a small amount 
[see Fig. 6(c)]. By varying this mass and monitoring the re- 
sulting change in resonance frequency, we produce the curve 
seen in Fig. 6(a). Extracting an effective mass ratio from the 
fit to this data gives a value of m^^ \/m — 0.3864. 

By examining Table Vll, we can see that the methods of 
Sader et al. and Cleveland et al. adapted for a doubly clamped 
beam provide very good approximations to the effective mass 
of the fundamental mode. This agreement suggests that mak- 
ing such an analogy is reasonable for this simple geometry. 



C. Effective mass of a string 

As the mode shape for a string is very simple, calculation 
of its effective mass ratio by inputting Eq. (39) into the EMI 
is straightforward and gives [28] 



Mode 




Label (n) 


EMI - FEM 


1 


0.4969 


2 


0.4966 


3 


0.4961 


4 


0.4955 


5 


0.4946 


6 


0.4933 


7 


0.4908 


8 


0.4911 


9 


0.4898 


10 


0.4883 



TABLE VIII: Values of m^f{ „/m calculated by solving the EMI using 
FEM for the first ten modes, n, of a simple string. 



^ - dx \u„[x)\ = -. (67) 

m L Jo 2 

Similar to the cantilever case, the effective mass is mode inde- 
pendent. The analytical result above can also be seen to agree 
well with FEM solutions of the EMI, which are summarized 
in Table Vlll. 

As with the doubly clamped beam, we can attempt to adapt 
the methods of Sader et al. and Cleveland et al. to determine 
approximations to the effective mass ratio of a string. 

Due to its geometry, it is very difficult to conceive of an an- 
alytical method by which one could calculate the spring con- 
stant of a string. However, it is possible to imagine applying 
the method of Sader et al. to a beam with an axial tension, 
which mimics a string at large tension values. We attempt 
to determine the effective mass of a string in this way, with 
the aid of FEM simulation. It should be noted that one must 
be very careful when applying such a force using FEM sim- 
ulation, as nonlinear effects must be considered in this case. 
Using this method, a value of mgff i/m — 0.4334 is found. 

Using the method of Cleveland et al., in which an added 
mass is applied at the center of the string, FEM simulation 
produces the data seen in Fig. 6. Fitting this curve with 
Eq. (62) gives a result of ;7Jeff.i //« = 0.4265. 

The values found using the methods of Sader et al. and 
Cleveland et al., while similar to each other, deviate signifi- 
cantly from the analytical value of 1 /2 we found in Eq. (67). 
This demonstrates a much greater discrepancy between exact 
and approximate values than what was observed when apply- 
ing similar methods to doubly clamped beams and cantilevers. 
This is due to the fact that the mode shape of the string is much 
more senstive to added mass, which is likely a result of the 
high tension in the device. This hypothesis is confirmed by 
observation of FEM mode shapes, as can be seen in Fig. 6(d), 
and by the fact that if we fit a segment of the curve which ex- 
cludes larger masses, we get a value that is closer to the expect 
effective mass ratio. 

It can be concluded that the methods of Sader et al. and 
Cleveland et al. provide an effective mass ratio for a string 
on the same order of that calculated using Eq. (61), however 
the approximate values deviate significantly from the exact 
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Method 




EMI - Analytical 


0.5000 


EMI - FEM 


0.4969 


Sader et al. - FEM 


0.4334 


Cleveland et al. - FEM 


0.4265 



TABLE IX: Values of meff,„/m for the fundamental mode of a simple 
string calculated using the different methods detailed in Sec. V. 



Mode Label 




{m,n) 


EMI - FEM 


(1.1) 


0.2498 


(2,1) 


0.2498 


(3,1) 


0.2498 


(1.2) 


0.2498 


(2.2) 


0.2498 


(4.1) 


0.2498 


(3,2) 


0.2498 


(1,3) 


0.2498 


(4,2) 


0.2498 


(5,1) 


0.2498 



TABLE X: Values of m^ff ,n„/m calculated by solving the EMI using 
FEM for the first ten modes of a rectangular membrane. The mode 
indices (m, n) match that of Fig. 4. 



ones. For this reason, the EMI should be used to calculate the 
effective mass of a string whenever possible, especially when 
precise results are impoitant. 



D. Effective masses of membranes 

1. Rectangularmembmnes 

We begin by calculating the effective mass ratio for a rect- 
angular membrane analytically by inputting Eq. (46) into the 
EMI giving 



Weff. 



m 



LxLy Jo Jo 



dxdy [\jf,n„{x,y)] 



1 

4' 



(68) 



This result is very simple and mode independent, as with the 
strings. In fact, we get the value that one would expect from 
multiplying the effective mass ratio for two strings, which is 
unsurprising as the mode shape found in Sec. IV was equiv- 
alent to multiplying the mode shapes of two strings. We also 
evaluate the EMI for a rectangular membrane using FEM sim- 
ulation, the results of which are summarized in Table X. It can 
be seen that these values agree very well with the analytically 
calculated value of 1/4. 



2. Circular membranes 

As with the rectangular membranes, we begin by calculat- 
ing the effective mass ratio of a circular membrane analyti- 
cally, inputting Eq. (52) into the EMI, which gives 



Mode Label 


'"eff.mn/™ 




(m, n) 


EMI - Analytical 


EMI - FEM 


(0,1) 


0.2695 


0.2693 


(1,1) 


0.2396 


0.2394 


(2,1) 


0.2437 


0.2435 


(0,2) 


0.1158 


0.1188 


(3,1) 


0.2357 


0.2356 


(1,2) 


0.1330 


0.1329 


(4,1) 


0.2254 


0.2253 


(2,2) 


0.1556 


0.1555 


(0,3) 


0.0737 


0.0753 


(5,1) 


0.2152 


0.2151 



TABLE XI: Values of meff.mn/™ calculated by solving the EMI both 
analytically and using FEM for the first ten modes of a circular mem- 
brane. The mode indices (m,n) match that of Fig. 5. 



Ka^ Jo 
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d<j) / sds [\j/,n„{s,(j))]\ (69) 
Jo 



This integral can be performed for both the cases m — and 
m > [80] and, using the boundary condition that Jm{oimn) = 
along with the Bessel function relation J,„-\{s) +Jm+\{'>) = 
2mJ„t{z) /z [63], we get the following effective mass ratios for 
circular membranes: 




if m = 0, 



[y,„+i («,„„)] ifm>0. 



(70) 



Again, the EMI can also be solved using FEM. The effec- 
tive mass ratio for a circular membrane calculated both ana- 
lytically and using FEM simulation are presented in Table XI. 
Upon comparing these values, it is evident that the analytical 
and FEM calculation match very closely. 

We note that, due to their poor results in calculating the 
effective mass ratio for a nanostring, the methods of Sader et 
al. and Cleveland et al. were not used to calculate the effective 
mass ratios for membranes, as the systems are analogous. 



VI. TORSIONAL RESONATORS 

Torsional nanomechanical resonators are very useful de- 
vices, as they are able to measure applied torques with excel- 
lent precision [12, 13, 26, 47-49]. By accurately calibrating 
these devices, studies of many physical phenomena, such as 
electric charge [5], magnetism [14, 18, 19], and superfluid- 
ity [25], can be explored. In this section, we will investigate 
how to properly calibrate such devices by using the procedure 
described above. 

Difficulties often arise when attempting to perform ther- 
momechanical calibration of a torsional resonator, as in most 
cases its motion cannot be described by a one-dimensional 
model. Asymmetries of the device [14] may call for a model 
which extends in a third dimension to describe the device's 
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(a) 
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■Wj 



w 



■W2 




FIG. 7: (a) Torsional device viewed from the negative z-axis. (b) De- 
flected mode shape for the fundamental torsional mode. Colorscale 
indicates displacement with blue as zero displacement and red as 
maximum displacement, (c) Fundamental mode of the torsional de- 
vice as viewed from the postive y-axis. The dashed line respresents 
the device's equilibrium position. 



geometry or may cause hybridization of torsional modes with 
other degrees of freedom, such as flexural modes [13]. How- 
ever, using the formalism described above, we can easily deal 
with such cases by the using the EMI, as it can be used to cal- 
culate the effective mass corresponding to any motion. This 
effective mass can then be used to properly calibrate the de- 
vice. 

It may seem awkward to attribute an effective mass (linear 
spring constant) to a torsional resonator, as we would expect 
that it would be described by an effective moment of inertia 
(torsional spring constant). However, what needs to be re- 
membered is that the definition of effective mass is arbitrary, 
and we have chosen in this paper to define it in such a way 
that it can be applied to any resonator. Therefore, it is possi- 
ble to thermally calibrate the motion of a torsional resonator 
using an effective mass. We must be careful here though, as 
the effective mass defined by the EMI coiTesponds to the dis- 
placement of the resonator and not its angular motion. Only 
once we have properly calibrated the device's displacement 
will we be able to deduce its angular motion. It is important 
to note that we could have alternatively chosen to calibrate the 
angular motion first and then use it to determine the device's 
displacement, but such a method differs from the calibration 
scheme described in this document. 

To demonstrate how the calibration method described in 
this paper can be used to calibrate the angular motion of a 
torsional device, we present the thermal calibration for the 
fundamental mode of the simple torsional resonator depicted 
in Fig. 7 as a case study. We will first calibrate its angular 
motion directly, then show how the exact same results can be 
achieved using the calibration method detailed above. 

Due to the symmetry of this device, only the angle of de- 
flection 0(f) must be known in order to fully describe the mo- 
tion of its fundamental mode, which is given by 



zi{x,t) =xtan0(f) Kx9{t), 



(71) 



where x and 0{t) are those defined in Fig. 7. Our utilization 
of the small angle approximation in Eq. (7 1 ) is justified by the 
fact that for thermal motion, the deflection angle is extremely 
small [13]. As our resonator exhibits oscillatory motion, we 
model 9{t) as a damped harmonic oscillator given by [81] 



e 



(72) 



where (Oi, Q, and /gff are the angular frequency, quality fac- 
tor, and effective moment of inertia for the fundamental mode 
of this device and T(f) is a time dependent torque applied to 
the system due to thermal fluctuations. We can determine the 
potential energy for the fundamental mode of the resonator: 



U = ^(ojpA jdx\ziix,t)\^ 

= ^(o^\e{t)\^pA jdxx^ 



(73) 
(74) 
(75) 



where we have defined the effective moment of inerta to be 
hn = pA jdxx^, which turns out to be normal definition of 
moment of inertia in this case. Note that in the above equa- 
tion we have made the approximation that the vast majority of 
the potential energy of the resonator will be due to the pad- 
dles of the torsion device and not the torsion rod itself. This 
approximation is valid for this device, as the moment of iner- 
tia of the paddles dominate that of the torsion rod. Applying 
the equipartition theorem to the above equation gives 



{U) = \hn(o^{e\t)) 



1 



-keT. 



(76) 



Upon inspection of Eqs. (72), (73), and (76), we see that 
is it identical to the corresponding equations in Sec. Ill if we 
take meff,„ ha, fit) x{t), a„{t) -> 0(f), and to,, ©i. 
Therefore, the thermal PSD for its angular motion is 



See{v) = 



2kH,sQ (v2-Vi2)2 + (VVl/e)' 



(77) 



or alternatively by defining an effective torsional spring con- 
stant as fCeff ~ Oi^/eff, we have [47] 



5ee(v) = 



2kBTv\ 



TZK^ffQ 



(V2 



' + {vvi/QY 



(78) 



We can then fit the experimentally obtained PSD data with a 
function of the form 



Svv{v) 



^VV ' 



pSeeiv). 



(79) 



This equation is identical to Eq. (17), except with a conver- 
sion factor j3, which allows for conversion of the theoretical 
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angular PSD (in units of rad^/Hz) to the experimental voltage 
PSD (units of V^/Hz). Once j3 is determined, we have all of 
the information we need to fully calibrate the angular motion 
of our data. 

It is easy to show that the calibration we just performed can 
be obtained by using the general method that was developed 
in the previous sections. We begin by determining Zi(x,t) for 
this resonator according to Eq. (3), which will be given by 



Zi{x,t) — ai{t)ui{x) — ai{t) 



2x 



w 



(80) 



Here w is the width of the torsional resonator as given in 
Fig. 7. Note that we have carefully chosen ui{x) such that 
|mi(x)| = 1 at its location of largest displacement located at 
X = zLw/2. The effective mass of this mode is then determined 
to be 



nigff 1 — pA / dx 



4x2 



'eff- 



(81) 



We can equate Eqs. (80) and (71), to obtain a relation be- 
tween 9{t) and ai (?) given by 0(f) = 2fli (f) /w. Using the re- 
lation between a device's motion and its PSD given by Eq. (8), 
we find that See{v) and S^^{v), which are the PSDs corre- 
sponding to 0(f) and a„(f), respectively, are related by 



Seeiv) ^ ^S,,{v) 



4 

vy2 



(v2-Vl2)2 + (vVi/e)' 



(82) 



2nH,fiQ (v2-v,2)V(vvi/e)' 



This gives us the exact same 5'ee(v) as we obtained calibrat- 
ing the angular motion directly. This tells us that we can 
equivalently describe the motion of this torsional resonator us- 
ing an effective mass or an effective moment of inertia. The 
only difference between the two is that we are measuring a 
different degree of freedom in each case, angular for the ef- 
fective moment of inertia and translational for the effective 
mass. 

We can also analyze our simple torsional resonator by ex- 
amining its effective mass ratio, as well as its effective mo- 
ment of inertia ratio /eff// where / = / dV p {r)r^ = pA f dxx^ 
is the conventional moment of inertia, by using methods sim- 
ilar to those used in Sec. V. Analytically calculating the ef- 
fective mass ratio for this torsional resonator by inputting 
Eq. (80) into Eq. (61) produces 



m 



1 /■«'/2 2 1 

- / dx |mi(x)| = -. 

W J-w/2 3 



(83) 



From this result, we can determine the effective moment of 
inertia for this system by using Eq. (81). Dividing this result 



by / = mn'2/12, which is the conventional moment of iner- 
tia for the rectangular paddle of the torsional resonator, gives 
hff/I = 1 for this device. This result is expected because, as 
mentioned before, in this situation the effective moment of in- 
tertia is defined to be identical to the conventional moment of 
inertia. Both the effective mass and moment of inertia ratios 
can also be calulated using FEM simulation, and are compared 
to the analytical values in Table Xll. 

We can also approximate the effective mass and moment 
of inertia ratios of our resonator by adapting the methods of 
Sader et al. and Cleveland et al. to be used on a torsional de- 
vice. For the method of Sader et al., we apply a known torque 
To to the resonator and measure its angular displacement 0o, 
which allows us to determine the torsional spring constant via 
Hooke's Law, Tq — —kQq. We can then approximate the ef- 
fective moment of inertia using the relation 4ff — fc/oi^. To 
perform these calculations, a FEM static solver is used, in 
which a force is applied to the two edges of the torsion pad- 
dle, equal in magnitude and opposite in direction, providing 
a torque and an angular displacement. This allows us to de- 
termine K, which we combine with (0\ obtained from FEM 
eigenfrequency analysis to calculate /gff using this method. 
The effective mass corresponding to this method can also be 
determined using Eq. (81). The results for both the effective 
mass and moment of inertia ratios using the method of Sader 
et al. are summarized in Table Xll. 

The method of Cleveland et al. can also be adapted for use 
on torsional resonators. In the analogous process, we observe 
the device's shift in resonance frequency due to an added mo- 
ment of inertia 4dd, which is described by 



1 

2n 



(84) 



'eff - 



^add 



To determine this frequency shift, we perform FEM simula- 
tion on the device, in which mass is added to the resonators 
paddles, providing an increased moment of inertia. Simula- 
tion data for such a procedure is presented in Fig. 8, along 
with a fit to Eq. (84), which allows us to extract an effective 
moment of inertia ratio and, via Eq. (81), an effective mass 
ratio, which are given in Table Xll. 

By investigation of Table Xll, we can see that the meth- 
ods of Sader et al. and Cleveland et al. provide a very good 
approximation of the effective moment of inertia and mass. 
However, it should be mentioned that while these methods 
work well for simple torsional devices, such as the one pre- 
sented here, they quickly break down as we move to hy- 
bridized modes, which have both torsional and flexural oscil- 
lations. For this reason, it is best to use the EMI to calculate 
an effective mass associated with a torsional device, and use 
it to calibrate the motion of the system. 

This simple example demonstrates that it is possible to ther- 
mally calibrate the motion of a torsional resonator using an 
effective mass, provided that we keep track of which degree 
of freedom we are measuring. This is an important result as 
it reinforces the fact that the method described in the previous 
sections can be used for any arbitrary oscillatory motion, be it 
flexural, torsional, or a complex combination of both. 
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0.0 0.2 0.4 0.6 0.8 1.0 
Normalized Added Moment of Inertia 

{4dd//) 

FIG. 8: Plot of normalized frequency versus normalized added mo- 
ment of inertia for the fundamental mode of the torsional resonator 
depicted in Fig. 7. The frequency values are normalized by di- 
viding by the fundamental frequency of each device with /add = 0. 
The added moment of inertia is normalized by dividing by the con- 
ventional moment of inertia / of the device. Solid lines are fit to 
Eq. (84). The effective moment of inertia ratio extracted from the fit 
is /eff// = 0.9948, which corresponds to an effective mass ratio of 
'Weff,l/'w = 0.3316. Both of these results can be found in Table Xll. 



Method 


mei=f,i/ra 


hn/l 


EMI ~ Analytical 


0.3333 


1.0000 


EMI - FEM 


0.3314 


0.9942 


Sader et al. - FEM 


0.3310 


0.9929 


Cleveland et al. - FEM 


0.3316 


0.9948 



TABLE XII: Values of in^ff i/m and /eff// for the fundamental mode 
of the torsional resonator depicted in Fig. 7 calculated using the dif- 
ferent methods detailed in Sec. V. 



VII. CONCLUSION 

This document provides a general procedure by which the 
thermal motion of any nano/micro-mechanical resonator can 
be properly calibrated via the equipartition theorem. It was 
shown that for any given resonator, a generic spectral density 
function can be used to describe its motion in the frequency 
domain. The only modification existing between different 
structural geometries was manifested by a single parameter, 
effective mass, which is a mode shape dependent quantity. 
Therefore, the effective mass was investigated in great detail 
for a number of nano/micro-resonator geometries, including 
cantilevers, doubly clamped beams, strings, membranes and 
torsional devices. The end result was a universal method by 
which any nano/micro-mechanical resonator can be thermo- 
mechanically calibrated, regardless of the complexity of its 
motion or number of degrees of freedom. By applying this 
method, the motion of such devices can be very precisely 
determined, allowing for use of nano/micro-mechanical res- 
onators for a number of different measurement applications, 
as well as in the investigation of fundamental scientific phe- 
nomena, with confidence and accuracy. 
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